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1  V  •  , 

-We-  discuss-nthe  initial  value  problem  in  one- dimensional  linear  visco¬ 
elasticity  with  a  step- jump  in  the  initial  data.  If  the  memory  kernel  is 

/  '  /  / 

sufficiently  smooth  on  [0,f>),  the  solution  exhibits  discontinuities 
propagating  along  characteristics  and  a  (higher  order)  stationary 
discontinuity  at  the  position  of  the  original  step-jump.  For  a  singular 
memory  kernel,  the  propagating  waves  are  smoothed  in  a  manner  depending  on  the 
nature  of  the  singularity  in  the  kernel,  but  the  stationary  discontinuity 
remains.  We  also  discuss  the  effects  of  these  phenomena  on  the  regularity  of 
solutions  with  arbitrary  initial  data. 
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SIGNIFICANCE  AND  EXPLANATION 


The  possibility  of  discontinuities  in  solutions  to  equations  which  model 
motions  of  viscoelastic  media  has  received  a  lot  of  attention.  In  this  paper, 
we  study  linear  initial  value  problems  with  jump  discontinuities  in  the 
data.  The  solutions  of  such  problems  can  exhibit  two  types  of  singularities: 
(i)  those  which  propagate  along  characteristics,  and  (ii)  stationary 
discontinuities  at  the  positions  of  the  original  jumps  in  the  data. 

In  particular,  we  are  interested  in  models  with  singular  memory 
kernels.  There  are  indications,  both  theoretical  and  experimental,  that 
singular  kernels  are  a  realistic  possibility  for  certain  viscoelastic 
materials.  Such  kernels  have  a  smoothing  effect  on  the  propagating  waves. 

The  degree  of  smoothing  depends  on  the  nature  of  the  singularity  in  the 
kernel.  Some  examples  lead  to  C  -smoothing  even  though  the  wave  speed  is 
finite.  There  are  also  cases  in  which  the  degree  of  smoothing  increases  in  a 
manner  proportional  to  time.  Kernels  with  very  weak  singularities  produce 
only  an  infinitesimal  gain  in  regularity.  Although  singular  kernels  tend  to 
increase  the  regularity  of  propagating  waves,  they  do  not  have  any  smoothing 
effect  on  the  stationary  discontinuities.  We  also  discuss  the  implications  of 
these  results  for  problems  with  arbitrary  initial  data. 
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ON  WAVE  PROPAGATION  IN  LINEAR  VISCOELASTICITY 


W.  J.  Hrusa1'2  and  M.  Renardy1'3 


1.  Introduction 

In  this  paper/  we  study  linear  wave  propagation  in  a  one-dimensional  viscoelastic 
medium.  That  is>  we  study  the  equation 

(1.1)  utt(x,t)  -  buxx(x,t)  +  m(t-T )  (u^fX/t )  -  u^fx.TjJdT,  xeDcR,  t>0  , 

where  b  is  a  nonnegative  constant  and  m  is  a  positive/  monotone  decreasing  kernel. 

(For  fluids  b  »  0,  while  for  solida  b  >  0.)  We  are  interested  in  initial  value  problems 

with  D  ■  R  and  discontinuous  initial  data.  More  precisely,  we  shall  study  situations 

where  u  is  identically  zero  for  t  <  0,  i.e. 

(1.2)  u(x,t)  -  0,  x  e  R,  t  <  0  , 

but  there  is  a  spatial  jump  in  the  instantaneous  values  of  u  and  ut  at  time  t  ■  0. 
Specifically,  we  look  at  the  following  types  of  initial  conditions 
(1.3a)  u(x,0)  -  j  sgn  x,  ut(x,0>  ■  0,  x  e  R  , 

(1.3b)  u(x,0)  -  0,  ut(x,0)  «  ~  sgn  x,  x  e  R  . 

(More  general  initial  conditions  can  be  treated  by  superposition.) 

There  have  been  a  number  of  studies  of  the  so-called  Rayleigh  problem,  i.e.  (1.1)  on 
the  half  line  D  «  [0,")  with  u  =  0  for  t  <  0  and  a  step-jump  in  the  boundary  value 

of  u  at  t  »  0  [2,  4-7,  11,  17,  18,  20-22,  24].  If  the  kernel  m  is  sufficiently 

smooth  (including  the  point  0),  then  the  discontinuity  at  the  boundary  will  propagate 
with  constant  speed  into  the  interior,  and  its  amplitude  will  decay  exponentially  in 
time.  On  the  other  hand,  Renardy  [22]  has  shown  that  certain  singular  kernels  lead  to 
smoothing  of  the  discontinuity.  More  precisely,  he  considered  the  Rayleigh  problem  with  a 
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class  of  kernels  for  which  m(t-T)  behaves  like  a  negative  power  of  (t-T)  near  zero,  and 


L  to 

i  showed  that  that  solution  is  of  class  C  for  t  >  0.  For  cases  involving  an  integrable 

[  singularity  in  m,  C  -smoothing  coexists  with  finite  speed  of  propagation  [22] . 

There  have  been  several  studies  of  linear  integrodif ferential  equations  in  Banach 
spaces  which  show  that  singular  kernels  can  have  a  regularizing  effect  on  solutions.  See, 

I  for  example,  [15]  and  the  references  therein. 

Hannsgen  and  Wheeler  [14]  have  shown  that  the  evolution  operator  associated  with  (1.1) 
has  a  certain  compactness  property  if  and  only  if  m(0+)  =  “.  This  indicates  that  singular 
kernels  lead  to  smoothing  of  solutions,  hut  does  not  characterize  the  degree  of 
smoothing.  In  fact,  as  we  shall  see  from  examples,  the  degree  of  smoothing  cannot  be 
specified  without  some  definite  assumption  on  the  nature  of  the  singularity  in  m. 

A  new  feature  that  arises  in  the  initial  value  problem  with  interior  discontinuities 
in  the  data  is  the  emergence  of  stationary  singularities  in  u.  Initial  conditions  (1.3a) 
and  (1.3b)  lead  to  solutions  with  discontinuities  in  u^  across  x  =  0. 

(A  notable  exception  occurs  for  b  «  0,  m(t)  =  Ke  wt,  and  initial  conditions  (1.3a),  in 
which  case  u  is  smooth  across  x  =  0  for  t  >  0.)  The  possibility  of  such  stationary 
singularities  for  hyperbolic  equations  with  memory  was  pointed  out  by  Greenberg,  Hsiao,  and 
HacCamy  [12],  [13],  [19],  but  does  not  seem  to  have  received  much  attention  otherwise.  In 
addition  to  the  stationary  singularity,  there  will  be  waves  propagating  to  both  sides.  The 
propagating  waves  are  given  by  the  same  kinds  of  expressions  as  those  obtained  in  the 
Rayleigh  problem,  and  consequently  the  same  regularity  results  apply. 

If  m  is  smooth  on  [0,<»),  there  will  be  propagating  discontinuities.  For  an 
appropriate  class  of  kernels  with  power-type  singularities,  the  propagating  waves  will  be 

OB 

of  class  C  .  Here  we  also  study  kernels  with  logarithmic  singularities,  for  which  we  find 
that  the  regularity  of  the  propagating  waves  increases  in  a  manner  proportional  to  time. 
Sufficiently  weak  singularities  in  m  (such  as  log-log)  produce  only  an  "infinitesimal" 
gain  in  regularity.  It  is  interesting  to  note  that  although  singular  kernels  have  a 
smoothing  effect  on  propagating  waves,  they  do  not  smooth  out  the  stationary 


discontinuities 


By  superposition,  these  results  can  be  used  to  study  the  regularity  of  solutions  with 

arbitrary  initial  data,  say  of  class  Hs.  This  is  discussed  briefly  in  the  last  section. 

If  m  is  smooth  on  [0,*),  then  the  strongest  singularities  in  u  are  always  the 

propagating  ones.  Roughly  speaking,  the  solution  preserves  the  smoothness  of  thi  data, 

00 

If,  however,  the  kernel  has  a  power-type  singularity,  the  propagating  waves  will  be  C 
and  there  will  be  a  fixed  (but  finite)  gain  in  regularity  determined  by  the  stationary 
discontinuities.  The  case  of  a  logarithmic  singularity  in  m  is  particularly  interesting. 
The  regularity  of  u  is  determined  initially  by  the  propagating  waves  and  increases  in  a 
manner  proportional  to  time  until  the  propagating  waves  and  stationary  singularities  are  of 
the  "same  strength".  Then,  the  stationary  singularities  dominate  and  regularity  does  not 
improve  any  further. 

In  order  to  focus  ideas  and  avoid  complicated  technical  hypotheses,  we  restrict  our 
attention  to  specific  classes  of  singular  kernels  that  are  given  by  infinite  sums  of 
exponentials.  (Kernels  of  a  similar  structure  appear  in  some  rheological  models  which  have 
been  motivated  by  molecular  considerations  [9],  [23],  [25] . )  It  is  clear  that  our 
arguments  can  be  applied  to  more  general  classes  of  kernels  with  similar  types  of 
singularities. 

Throughout  this  paper,  the  operations  of  differentiation,  convolution,  and  Laplace 
transformation  should  be  interpreted  in  the  sense  of  tempered  distributions.  (See,  for 
example,  [10].)  We  shall  not  dwell  on  this  point.  However,  we  note  that  the  classical  and 
distributional  definitions  agree  for  sufficiently  regular  functions.  Moreover,  since  the 
distributional  definition  of  an  operation  typically  involves  performing  the  operation  on 
smooth  test  functions,  there  is  generally  no  difficulty  with  changing  the  order  of 
"limiting"  operations.  In  particular,  the  order  of  "differentiation  with  respect  to  x" 
and  "Laplace  transformation  with  respect  to  t"  can  always  be  interchanged  in  the 
distributional  setting. 

Acknowledgement i  We  are  indebted  to  YUriko  Renardy  for  some  valuable  suggestions. 
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2.  Initial  Value  Problems  with  Step-Jump  Initial  Conditions 

In  this  section  we  investigate  solutions  of  (1.1)  with  D  “  R.  Me  assume  throughout 
that  m  is  nonnegative  and  nonincreasing,  but  m  ?  0,  and  that  b  >  0.  We  are  interested 
in  situations  where  u  £  0  for  t  <  0,  but  initial  conditions  (1.3a)  or  (1.3b)  hold  as 
t  •*  0+.  To  solve  (1.1),  we  introduce  Laplace  transforms  with  respect  to  the  time  variable: 

(2.1)  u(x,X)  =  e  ^tu(x,t)dt 
Observe  that  (1.3a),  (1.3b)  can  be  subsumed  under 

(2.2)  u(x,0)  =  j  a  sgn  x,  ufc(x,0)  -  ~  0  sgn  x  . 

It  is  easy  to  see  that  (1.1),  (1.2),  (2.2)  lead  to  the  equation 

(2.3)  X2  u(x,X)  -  aX  sgn  x  -  B  sgn  x  =  ( b+m(0  )-m(A )  )uxj((x,X  ) 

*  A  ,ce  f 

(We  remark  that  "m (0)  -  m(A>"  should  be  interpreted  as  m(t)(1-e  )dt,  which  may 

converge  even  if  m  has  a  nonintegrable  singularity  at  0,  i.e.  m(0)  and  m(A )  do  not 
make  sense  separately.) 

For  fixed  X  this  is  an  ordinary  differential  equation  which  is  easily  solved.  We 
find  that 

(2.4)  u(x,X)  =*  (aX+S){— 1-  sgn  x - ~  H(x)exp( -Ax//b  +  m(0>  -  m(X)  ) 

2A 1  2\ 


+  — rr  H(-x)exp( Ax//b  +  m(0)  -  m(X)  )} 

2XZ 

where  H  denotes  the  Heaviside  step  function.  The  solution  of  (1.1),  (1.2),  (2.2)  is 
obtained  by  taking  the  inverse  Laplace  transform  of  u.  Although  this  procedure  has  been 
purely  formal,  it  can  be  justified  rigorously  under  rather  mild  assumptions  on  m.  See 
[20]  for  a  discussion  of  this  in  a  similar  situation. 

For  x,  t  >  0  we  have 

(2.5)  u ( x , t )  =  j  (a+$t)  -  /***“  +  ^~)exp(-Xx//b  +  m(0)  -  m(X)  ) 

where  y  is  any  positive  number,  and  the  integrand  is  made  single-valued  in  the  usual 

way.  A  similar  expression  holds  for  x  <  0.  In  fact,  u(-x,t)  ■  -u(x,t)  for  all 

*  1 
x,  t  >  0.  Observe  that  if  m(0)  -  m(A)  is  well-defined  (in  particular,  if  m  e  L  ( 0 ) 


or  if  tm(t)dt  <  “  and  /*  m(t)dt  <  * ) ,  than 
(i)  Re(m(0)  -  m(A))  >  0  for  Re  A  >  0  > 

A  A 

(11)  lm(m(0)  -  m(A))  >  0  for  Re  A  >  0,  In  X  >  0  t 

A  A 

(111)  lm(ra(0)  -  m(A))  <  0  for  ReX>0,  ImA<0  , 

by  virtue  of  our  sign  conditions  on  n. 

He  note  that  if  m  belongs  to  l/(0,<*)  then  the  solution  of  (1.1),  (1.2),  (2.2)  has 
finite  speed  of  propagation.  Indeed,  for  m  e  I^tO,"),  the  Riemann-Lebesgue  lemma  Implies 
that  m(A)  ♦  0  as  1  M  1  >  0.  Consequently,  the  contour  of  integration  in  (2.5)  can 

be  closed  by  a  circle  on  the  right  if  x  >  t  •  /b  ♦  m(0 ) .  Since  the  integrand  has  no 
singularities  with  Re  X  >  0,  it  follows  from  Cauchy's  theorem  that  the  integral  in  (2.5) 


vanishes  for  x  >  t  •  /b  +  m(0) .  (See  (41.)  The  situation  is  analogous  for 
x  <  -t  •  t/b  +  m(  0 ) . 


A.  Regular  Kernels 

Our  objective  here  is  to  describe  the  regularity  of  solutions  when  the  kernel  ra  is 
smooth  on  (0,“).  In  order  to  highlight  the  main  ideas  and  avoid  repeated  hyptheses,  we 
assume  that 

*  (k)  1 

(2.6)  m  e  C  (0,-),  m'  '  e  L  (0,-)  V  k  •  0,1,...  , 

i.e.  that  m  and  its  derivatives  of  all  orders  are  continuous  on  [0,0")  and  integrable. 

*  1 
It  follows  that  m(X)  can  be  approximated  asymptotically  by  a  power  series  in  —  as 

X  ♦  «,  Re  X  >  0.  This  can  be  used  to  study  smoothness  properties  of  the  solution.  (See 

C6J.) 

We  first  discuss  the  situation  for  x  >  0.  Away  from  the  line  x  •  t  •  /b  +  in(0 ) , 

u  is  of  class  C  .  Across  this  line,  u  sustains  a  jump  discontinuity  in  case  (1.3a), 
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while  in  caee  (1.3b),  u  is  continuous  but  ux  and  ^  jump.  The  amplitudes  of  these 


jumps  decay  exponentially  in  time.  In  fact,  for  x  >  0  the  solution  of  (1.1),  (1.2), 

(2.2)  can  be  written  in  the  form 

(2.7)  u ( x , t )  -  j  (a  +  8t) 

,  -  ,/,  _ 2 - 

-■5- a  exp[-m(0  )x/2(b  +  ■<0)r/  J  •  H(.t  -  . - r - J 

2  /b  +  m(0> 

-  |  [8 - ^J0)ay_ - ?gm(0_) fx  jexp[-m(0) x/2/b  +  m(0»  • 

2  2(b+m(0))J/  8(b  +  m(0) )a/z 

X  X 

(t  -  , - „ - )  •  H(t  -  , - x - ■) 

/b  +  m(0 )  rb  +  m(0 ) 

+  F(x,t)  , 

where  F  is  of  class  C1  across  the  line  x  *  t  •  h  +  m(0)  and  of  class  C*  elsewhere, 
and  F(x,t)  *  0  for  x>t*^+m(0).  Expressions  such  as  (2.7)  are,  of  course,  well 
known.  (See,  for  example,  [2).) 

The  situation  is  completely  analogous  for  x  <  0.  Thus,  the  behavior  of  u  away  from 
the  line  x  «>  0  is  quite  similar  to  the  case  of  a  linear  wave  equation  with  frictional 
damping.  However,  across  x  *  0  solutions  of  (1.1),  (1.2),  (2.2)  generally  exhibit 
discontinuities  in  uxx  which  do  not  occur  for  equations  without  memory.  (An  exception 
arises  in  the  special  case  m(t)  =  Fe  b  *  0  and  initial  conditions  (1.3a)  hold.) 

Discontinuities  of  this  type  have  been  discussed  in  a  series  of  papers  (concerning  Riemann 
problems  for  first  order  hyperbolic  equations  with  memory)  by  Greenberg,  Hsiao,  and  MacCamy 
[12],  (13),  (19),  but  do  not  seem  to  have  received  much  attention  otherwise. 

/  *  00 

Away  from  the  lines  x  »  0,  x  ■  ±  t  •  rb  +  m(0),  u  is  of  class  C  .  Moreover,  for 
t  >  0,  all  derivatives  of  u  have  (finite)  right  and  left  limits  as  x  ♦  0.  It  is  not 
difficult  to  show  that  u,  ut,  ux,  utt,  utx  (and  uxt)  are  continuous  across  x  =  0.  To 
investigate  the  possibility  of  a  stationary  discontinuity  in  u^,  we  set 

(2.8)  J(t)  »  u^tO^t)  -  uxx ( 0™ , t )  V  t  >  0  . 

It  follows  from  (2.4)  and  a  straightforward  computation  that 


(2.9) 


J(X) 


am(O)m(X)+tt(b+m(O))(lm(X)-m(O))+0(b+m(O))m(X) 

(b+m(0))2(b+m(0)-m(X>) 


We  note  that  the  distributional  derivative  uxx  is  the  sum  of  a  function  and  a 
distribution  supported  on  the  lines  x  -  ±  t  I b  +  m(0 )  (which  arises  from  the  propagating 
jumps  described  above).  The  part  of  uxx  which  is  supported  on  these  lines  makes  no 
contribution  to  J(t)  and  has  been  "subtracted  off"  in  our  derivation  of  (2.9).  If  one 
were  to  drive  an  expression  for  J(X)  by  simply  taking  right  and  left  derivatives  as 
x  *  0  in  (2.4),  then  this  expression  would  contain  additional  terms  corresponding  to  a 
distribution  with  support  at  0. 

The  first  question  that  we  address  regarding  J  is  whether  or  not  there  are  any 

(nontrivial)  situations  in  which  J  vanishes  identically.  It  is  clear  from  (2.9)  that  if 
2  2 

a  +  0  f  0  and  J  5  0,  then  m  must  be  an  exponential  of  the  form 
(2.10)  m(t)  -  Ke_Wt  ,  K,  W  >  0  . 

(Recall  our  sign  assumptions  on  m.)  Further  examination  of  (2.9)  reveals  that  under 
initial  conditions  (1.3a),  J  3  0  if  and  only  if  b  «  0  and  m  is  of  the  form  (2.10). 
Moreover,  for  initial  conditions  (1.3b)  there  are  no  cases  in  which  J  vanishes 
identically.  If  values  of  a,  0  other  than  those  corresponding  to  (1.3a),  (1.3b)  are 
considered,  then  there  will  be  additional  situations  in  which  J  S  0,  eg.  a  »  0  =  1, 
b  -  j.  K  ■  1,  u  ■>  2.  (These  can  easily  be  characterized  using  (2.13),  (2.14),  (2.15) 
below. ) 

We  note  that  JSC  (0,“).  An  expression  for  J(0)  can  be  found  by  computing 

lira  X J ( X ) .  The  outcome  of  this  simple  computation  is 
X»- 


(2.11) 


J(0) 


cm(0 )  +a(b+m(0))m' (0 )+0( b+m(0 ) )m(0 ) 
(b+m(0)  )3 


Observe  that  if  b  «  0  and  initial  conditions  (1.3b)  hold,  then  J(X)  has  a  pole  (or 
worse)  at  X  »  0  and  consequently  one  should  not  expect  decay  of  J(t)  as  t  +  •».  On  the 

other  hand,  our  sign  conditions  on  m  imply  that  if  b  >  0  or  if  b  =  0, 

<  * 
tm(t)  8  L  (0,“),  and  initial  conditions  (1.3a)  hold,  then  J(X)  has  no  poles  with 

Re  X  >  0,  and  one  does  expect  that  J(t)  +0  as  t  *  «. 
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By  virtue  of  (2.6),  our  sign  conditions  on  m,  and  the  Riemann-Lebesgue  lemma,  it 

follows  from  (2.9)  and  Proposition  2.3  of  [161  that  J  S  L^(0,“)  if  b  >  0  or  if 

b  =  6  “  0  and  t m(t>  e  b'fO,-).  By  the  same  kind  of  reasoning,  we  find  that  J'  also 

belongs  to  l/(0,*>)  under  either  of  the  above  conditions.  (In  fact,  so  do  all  higher 

order  derivatives  of  J. )  We  therefore  conclude  that  if  b  >  0  or  if  b  ■  6  “  0  and 

tm(t)  e  L^O,®),  then  J(t)  *0  as  t  ♦  ». 

A  similar  argument  can  be  used  to  handle  the  case  when  b  =  0  and  initial  conditions 

*  .  —  1  , 00  —  1 

(1.3b)  hold.  Examination  of  (2.8)  suggests  that  J(t)  ♦  [(m)’(0)]  »  Tm(T)dt)  as 

t  ♦  “.  The  idea,  then,  is  to  show  that  J»,  J1  6  l'(0,"“),  where 
(2.12)  J*(t)  =  J(t)  -  [ (m) ' (0 ) )-1  V  t  >  0  . 

If  we  assume  that  b  =  0,  tm(t)  e  L  (0,®),  and  initial  conditions  (1.3b)  hold,  then  it 
follows  from  Proposition  2,3  of  [16J  and  some  routine  computations  that  J*,  e  L*(0,»), 

fOO  mm  \ 

and  consequently  J«(t)  ■*  0  as  t  ♦  ®,  i.e.  J(t)  +  “wg  T»(t)dT)  as  t  •»  ®. 

The  asymptotic  behavior  of  J  can  thus  be  summarized  as  follows.  If  b  >  0,  then 
J(t)  +  0  as  t  •»  «*.  If  b  =  0  and  tm(ti  e  L^O,®),  then  J(t)  ♦  -B(Jq  Tm(T)dt)-1  as 
t  ■»  “.  This  last  result  suggests  that  if  b  ■  0  and  m(  t  >  ~  t  * 1 +0  *  as  t  ♦  "  with 
0  <  6  <  1,  then  J(t)  *  0  as  t  ♦  ®>  however,  we  have  been  unable  to  prove  it.  See  the 
paper  of  Chu  [6]  for  an  interesting  discussion  of  the  asymptotic  behavior  of  u  as  t  ♦  ® 
along  lines  of  the  form  x  =  a*t,  a  >  0* 

Since  the  Laplace  transform  of  a  function  (or,  more  generally,  a  distribution)  with 
compact  support  is  necessarily  an  entire  function,  if  follows  immediately  from  (2.9)  that 
if  b  *  0  and  3^0  then  the  discontinuity  in  uxx  persists  for  all  t  >  0.  This  is 

also  the  case  in  a  number  of  other  situations.  In  particular,  we  have  been  able  to  show 

2  2 

that  the  support  of  J  is  unbounded  if  b>0,  a3  *  0,  and  a  +3  f  0,  or  if  m  has 
2  2 

compact  support  and  a  +6  ^0.  The  proof  is  straightforward,  but  not  entirely 

trivial .  We  omit  the  details.  (See  additional  comments  concerning  completely  monotonic 
kernels  in  part  E  below.) 
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An  interesting  expression  for  J  can  be  obtained  by  expanding  lb  +  m(0)  -  m(X)] 
in  powers  of  m(X)/(b  +  m(0))  and  formally  inverting  (2.9)  as  a  sum  of  iterated 
convolutions.  The  series  for  J  obtained  in  this  manner  converges  in  L 1 ( 0 )  if  b  >  0 
and  in  [0 ,“)  if  b  -  0.  This  provides  an  alternate  proof  of  the  integrability  of 

J  when  b  >  0 . 

The  situation  is  particularly  simple  for  a  single  exponential  of  the  form  (2.10).  As 
noted  above/  for  b  ■  0  and  initial  conditions  (1.3a)  we  have  J  =  0.  Moreover,  in  this 
case  it  is  straightforward  to  show  that  higher  order  derivatives  of  u  are  continuous 
across  x  -  0.  For  b  >  0  and  initial  conditions  (1.3a),  we  conclude  from  (2.8)  that 

(2.13)  J< t )  «  _MH__  expt-ybt/(b+r)l  , 

(btr  r 

where 

(2.14)  r  »  m(0)  «  K/y  . 

A  similar  expression  holds  for  initial  conditions  (1.3b),  namely 

(2.15)  J(t)  -  — exp[-ubt/(b+r  )  1 

(b+r ) 

It  is  interesting  to  observe  that  for  b  »  0  and  the  kernel  (2.10),  the  initial  value 
problem  (1.1),  (1.2),  (2.2)  is  equivalent  to 

(2.16)  utt<x,t)  +  uut(x,t)  =  ^  u^tx.t)  +  wB  sgn  x  , 

with  initial  conditions  (2.2).  We  see  from  equation  (2.16)  that  if  &  t  0  then  the 
solution  must  have  a  singularity  across  x  «  0.  Recalling  that  ut  and  utt  are 
continuous  across  x  =  0  we  can  deduce  directly  from  (2.16)  that  J(t)  “  -y^6/X,  which 
agrees  with  our  previous  expressions  for  this  case. 

8.  Kernels  with  Power-Type  Singularities: 

In  [22],  Renardy  studied  the  Rayleigh  problem  for  the  class  of  singular  kernels 

°°  P 

(2.17)  m(t)  *£ent  ,p>—  . 

n»1 

As  t  ♦  0+,  m(t)  behaves  like  t"1/,p.  Accordingly,  m  is  integrable  if  p  >  1.  For 
p  >  1,  we  have 
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(2.18) 


1 


m(X)  -  l  - 

n»1  X+np 


while  for  -j  <  p  <  1  we  interpret  m(0)  -  m(X)  as 


(2.19)  /“<1-e"Xt)m(t)dt  -  7  - - -  . 

°  n-1  np(Xtnp) 

Using  the  results  of  (221,  we  can  immediately  draw  the  following  conclusions  regarding 

(1.1),  (1.2),  (2.2)  with  m  given  by  (2.17):  For  p  >  1  the  solution  u  is  analytic 

/ - - - 

away  from  the  lines  x  ■  0,  x  “  ±  t  •  »b  +  m(0),  it  is  of  class  C  across  the  lines 
x  ■  ±  t  •  ^b~+~m<o7 ,  and  it  propagates  with  finite  speed;  for  j  <  p  <  1,  u  is  analytic 
away  from  the  line  x  -  0.  Thus  the  singularity  in  m  smooths  out  the  propagating 
waves.  However,  it  does  not  smooth  out  the  stationary  discontinuity  in  u^. 

By  deforming  the  contour  of  integration  in  (2.5)  as  indicated  below,  it  is  straight¬ 
forward  to  show  that  u  and  its  derivatives  of  all  order  have  (finite)  limits  as  x  *  0+ 
for  t  >  <3. 


la  X 


Fig.  1.  Contour  of  Integration 

For  p  >  1,  this  contour  can  be  used  provided  that  0  <  x  <  t  •  /b  +  in(0 ) ;  if 
j  <  p  <  1,  it  can  be  used  for  all  x,  t  >  0.  The  situation  is  similar  as  x  *  0_. 
Moreover,  u,  ut,  ux,  u^,  and  \x^x  are  continuous  across  x  -  0  for  t  >  0. 
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Our  previous  expression  for  J(X)  is  not  valid  under  the  present  circumstances. 


Nevertheless/  it  is  not  too  difficult  to  show  that  uxx  is  discontinuous  across  x  *  0. 

The  argument  proceeds  roughly  as  follows.  Using  (2.4)  we  find  that 

(2.20)  u  (x,X)  ♦  — -  as  x  0+  . 

2<b+m(0)-m(X)) 

(We  note  that  there  is  no  problem  with  interchanging  the  order  of  differentiation  with 
respect  to  x  and  taking  Laplace  transforms  -  even  in  the  classical  sense  -  for  x  >  0.) 
Suppose  that  uxx  is  continuous  across  x  •  0  for  t  >  0.  Then,  it  is  easy  to  see  that 
uxx(0,t)  “  0  for  all  t  >  0  and  consequently  uxx(x,t)  *  0  uniformly  for  t  in  compact 
subsets  of  (0,*)  as  x  ♦  0+.  It  follows  that  the  expression  on  the  right-hand  side  of 

(2.20)  must  be  the  Laplace  transform  of  a  distribution  with  support  at  0,  i.e.  it  must  be 
a  polynomial  in  1.  However,  it  is  easy  to  show  that  the  expression  in  question  is  not  a 
polynomial,  and  we  conclude  that  u^  is  discontinuous  across  x  -  0  (provided,  of 
course,  that  a2  ♦  8  2  +  0.) 

The  discontinuity  in  u^  persists  for  all  t  >  0.  Indeed,  using  the  contour  of 
Figure  1  in  the  inversion  integral,  one  readily  sees  that  the  jump  J(t)  (defined  by 
(2.8))  is  analytic  in  t  >  0.  Therefore,  since  3  t  0,  it  cannot  vanish  identically  on  an 
Interval  of  the  form  (T,«). 

One  can  also  argue  as  follows  to  show  that  the  support  of  3  is  unbounded:  If  J 
had  compact  support. than  the  expression  on  the  right-hand  side  of  (2.20)  would  be  the 
Laplace  transform  of  a  distribution  with  compact  support  and  hence  an  entire  function  of 
X.  Examination  of  (2.18)  and  (2.19)  reveal  that  this  cannot  be  the  case  (unless,  of 
course ,  a  *  U  *  0 ) . 

C.  Kernels  with  Logarithmic  Singularities: 

An  example  of  kernel  which  behaves  nicely  at  infinity  and  has  a  logarithmic 
singularity  at  zero  is  given  by 

“  _  i> 

(2.21)  m(t)  -  l  e~e 

n»0 

To  verify  this,  we  note  that  the  asymptotic  behavior  as  t  0+  does  not  change  if  the  sum 


is  approximated  by  an  integral.  More  precisely,  using  the  estimate 


(2.22) 


f(v)dv  -  l  f (n) |  <  ]*|f • (v) |dv 


n-1 


with  f(u)  =  exp(-e  t),  we  find  that 
(2.23) 


We  then  observe  that 
(2.24) 


!/;  eVtdv  -  l  e'6”^  <  e'6  Vt  >  0  . 

n*  1 


v 

.»  -e  t .  rm 

In  e  dv  *  L  e 


dv 


/Int  e"*  dv 


/«  — e  fO  — « 

“  *0  *  dV  +  lint  *  dv  • 

The  last  term  in  (2.24)  clearly  has  a  logarithmic  singularity  as  t  ♦  0+.  We  note  also 
that  m  e  ( 0,“). 

It  follows  from  (2.21)  that 


(2.25) 


m(A)  -  l  . 


n«1  A+e 


To  determine  the  asymptotic  behavior  of  m(A)  as  |A|  ♦  •>,  we  set  y»ReA,4i»Imi, 
and  observe  that 
(2.26) 


ip  ~~TT  ”  T  4n(1+A, 


'°AteV  * 

Using  (2.26)  and  the  estimate  (2.22)  with  f(v)  *  (A+«V)  \  we  deduce  that 
(2.27)  m(A)  «  j  in(1+A)  +  0(|if>l  ')  as  !<i|  ♦  *•  . 

Therefore,  we  have 

(2.2B)  -  ■■  ■  ------ - * - - T7S-  +  - —  +  0(1) 


£ 7 


m(0  )  -  m(A ) 


(b  +  m(0))1/2  2(b  +  m(0))3/2 


and  consequently 

(2.29)  exp[-Ax//b  +  m<0)  -  m(A)) 


exp[ - ^ — _].expl — ^i^illiL—J.exptOd)] 


1/2 

(b  +  m(0) )  2(b  +  m(0 ) ) 

uniformly  for  x  in  bounded  sets  as  |v|  ♦  -. 
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Since  m  belongs  to  L 1  (0,“) ,  the  integral  in  (2.5)  vanishes  if  x  >  t  *  /b  *  m(0). 


For  0  <  x  <  t  •  »b  +  m(0 ) ,  we  can  deform  the  contour  of  Integration  as  shown  in  Figure  1. 
Along  this  contour,  the  integral  is  absolutely  convergent  and  remains  absolutely  convergent 
if  x  and  t  are  extended  into  the  complex  plane.  We  can  therefore  easily  conclude 


that  u  is  analytic  in  the  region  0<x<t*»bem(0). 

The  regularity  of  u  across  the  line  x  "  t  •  /b  +  m(0)  is  directly  related  to  the 
decay  of  expt-x  4n(1+X)/2(b  +  m(0))3/^2]  as  jX|  ♦  ■.  Since  this  term  behaves  like  a 
negative  power  of  X  with  exponent  proportional  to  x,  we  can  deduce  from  (2.5),  (2.27), 


and  (2.29)  that  the  regularity  of  u  across  the  line  x  “  t  •  »b  +  m(0 )  increases  in  a 
manner  proportional  to  x  (or  equivalently  proportional  to  t). 

It  is  interesting  to  observe  that  discontinuities  in  the  initial  data  become 
continuous  instantaneously.  More  precisely,  the  solution  of  (1.1),  (1.2),  (1.3a)  is 


continuous  across  x  »  t  •  rb  +  m(0)  for  t  >  0,  and  the  solution  of  (1.1),  (1.2),  (1.3b) 
is  of  class  across  x  «  t  •  i/b  +  m(0 )  for  t  >  0.  This  follows  immediately  from  the 

decay  properties  of  the  integrand  in  (2.5). 

The  situation  is,  of  course,  completely  analogous  for  x  <  0.  Arguing  as  before,  one 
can  show  that  u  and  its  derivatives  of  all  orders  have  (finite)  right  and  left  limits  as 
x  -»  0  for  t  >  0,  Moreover  u,  u^,  u„,  u^j.,  and  u^x  are  continuous  across  the  line 

x  »  0  for  t  >  0.  However,  uKX  sustains  a  Jump  discontinuity  across  this  llnej  as 

2  2 

before  the  discontinuity  persists  for  all  t  >  0  if  o  +B  ¥  0. 


D.  Kernels  with  Weaker  Singularities 

As  the  previous  examples  show,  the  regularity  of  u  across  the  lines 


x  «  t  t  •  /b  m(0)  depends  crucially  on  the  behavior  of  m(t)  as  t  ♦  0+.  Weak 
singularities  in  m  generally  lead  to  less  smoothing  than  strong  singularities  do.  In 
fact,  a  sufficiently  weak  singularity  in  m  will  produce  a  gain  in  regularity  which  is 

effectively  Infinitesimal,  we  illustrate  this  with  the  kernel 

•  n 

(2.30)  m(t)  »  £  exp(-t*e®  ) 

n»1 


which  has  a  log- log  singularity  at  0 


As  before,  the  behavior  of  ra  near  0  and  the  behavior  of  m  near  infinity  do  not 

change  if  the  sum  is  approximated  by  an  integral.  It)  verify  the  nature  of  the  singularity 

in  m,  we  observe  that  for  t  <  1# 

v 

I  expt-te®  )dv 


(2.31) 


r 

Jo 


in  in  t  ,  (e  +in  t),. 

1  expt-e  ]dv 


,  r"  r  (e  tin  t),. 

+  ^tn  jin  t|exp!_e  )dV  ‘ 


The  integrand  in  the  first  term  on  the  right-hand  side  is  bounded  above  and  below  by 
positive  constants,  hence  this  term  behaves  like  in | In  t|  as  t  ♦  o+.  With  the  substi¬ 
tution  n  3  v  -  in | in  t|,  the  second  Integral  becomes 

, „  ...  ,*•  ,  (e'1-1)»|in  t|,, 

(2.32)  J0  expl-e  ')dn 

which  is  clearly  bounded  as  t  ♦  0+.  Using  (2.22),  it  is  easy  to  see  that  the  difference 
between  m(t)  and  the  integral  on  the  left-hand  side  of  (2.31)  remains  bounded  as  t  ♦  0+. 
We  have  thus  verified  that  m  has  a  log-log  singularity  at  0. 

It  follows  from  (2.30)  that 


(2.33) 


m(X)  »  l 


n»1  X+exp(e  ) 

and  again  the  asymptotic  behavior  as  |X|  ♦  •  does  not  change  if  the  sum  is  approximated 
by  an  integral.  (As  in  part  C,  the  error  in  this  approximation  is  0(|ii|  1 ) . )  For  real 
X  let  us  set 

dv 


(2,34) 


X+exp(e  ) 


With  the  substitution  w  *  exp(e  ),  we  find  that 

dw 


(2.35) 


ItX>  ”  K  <X+w)w  in 


/ 


dw 


dw 


e  (X+w)w  in  w  ^X  (X+w)w  in  w 

Expanding  (X+w)  1  in  powers  of  w/1  and  X/w,  respectively,  and  then  making  the 
substitution  v  »  in  w,  we  see  that 
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(2.36) 


i(X)-|  l  <-i>Vn  /x  5! — dw 

X  n-0  e  tn  w 


l  (-1  )nXn  /' 


n.  n  f*  dw 


n-0 


X  n+2, 
w  tn  w 


»  j  tntn  X 


+  1  I  <-i)nx'n  j*n  x  s!ldv 

X  ,  '1  v 

n-  i 


+  1  v  <-nnin  r  ay 

X  *•  ^in  X  (n+1)v 

n-0  ve 


■1  tntn  X  +  l  ■  [El (n  in  X)  -  El(n)] 
n”  1  Xn 

at 

+  I  (-1)"xnEi((n+1)tn  X)  , 
n-1 


where  El  denote*  the  exponential  Integral.  (See,  for  example,  page  228  of 
eecond  tent  in  the  laat  line  la  0<-^-)  aa  X  *  »,  and  UBlng  formula  5.1.51 


[1]  we  find  that 
(2.37) 


l  (-1)nXnEl((n+1)tn  X)  -  *  o( - - - -J  . 

n-1  a  tn  A  XUn  x,z 


He  therefore  conclude  that 


(2.38) 


i 


dv 


0 - V - X  ln*n  X  *  XT nT  * 

0  X+exp(e  )  A  tn  A  X(tn  X) 


as  X  ♦  •.  By  analytic  continuation,  (2.38)  holds  for  complex  X  in  a  right 
Since  the  difference  between  m(X)  and  the  Integral  in  (2.38)  is  0<  | ip |  '), 

(2.39)  m(X)  -  j  tntn  X  +  0(|*|_1) 

aa  It!  *  •  for  suitable  y  >  0.  (Recall  that  t  ■  It*  X  and  p  -  Im  X.) 

Away  from  the  lines  x  «  ±  t  *  rb  +  m(0),  the  smoothness  properties  of 
exactly  the  same  as  in  part  C.  It  follows  from  (2.39)  that  the  integrand  in 
written  in  the  form 

x  . 

(2.40)  exp[x(t  -  - - x - )J  •  G(x,X>  •  [y  + 

/b  +  m(0)  X* 


(11.)  The 
on  page  231  of 


half  plane, 
we  have 

u  are 
(2.5)  can  be 
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where  G(x,A )  decays  like  a  negative  power  of  fin  A |  with  exponent  proportional  to  x. 

The  factor  G(x,A)  leads  to  an  "infinitesimal*  gain  in  regularity  acroea  the  line  x  • 

1 - * - 

t  •  »b  +  m(0).  The  situation  is,  of  course,  completely  analogous  for  x  <  0. 

It  is  interesting  to  observe  that  for  t  sufficiently  large,  the  solution  of  (1.1), 
(1.2),  (1.3a)  is  continuous  across  the  lines  x  »  t  t  •  /b  +  m(0 )  and  the  solution  of 
(1.1),  (1.2),  (1.3b)  is  of  class  across  these  lines.  (Indeed,  the  appropriate 

inversion  integrals  converge  absolutely  if  fxf  is  sufficiently  large.)  He  conjecture 
that  the  above  continuity  properties  actually  hold  for  all  t  >  0,  but  have  been  unable  to 
prove  it.  The  difficulty  lies  in  obtaining  a  suitable  error  estimate  for  our  approximation 
of  m(A).  One  can  show  that  the  principal  term  in  (2.39)  leads  to  a  function  which  is 
continuous  across  x  ■  t  •  Jh  +  m(0 ) .  However,  even  though  the  remainder  decays  more 
rapidly  than  the  principal  term  as  { X  {  ♦  “,  its  decay  rate  is  not  sufficient  to  guarantee 
the  desired  continuity.  Thus,  more  detailed  information  is  needed. 


We  can  prove  contuity  of  the  wave  front  for  the  kernel  given  by 


(2.41) 


m(t)  ■  /”  exp(-t  *  ee  )dv  , 


which  also  has  a  log-log  singularity  at  0.  In  this  case,  we  have 


(2.42) 

and  consequently,  by  (2.37) 

(2.43) 


m(A)  «  / 


dA 


0  v 

X  *  exp(e  ) 


;<x> "  x tntn  A  +  At  ♦  °(- 


A(in  X)- 

as  | gi  |  ♦  ".  Once  again,  it  follows  that  there  is  only  an  infinitesimal  gain  in  regularity 

_____ 

across  the  lines  x  »  ±  t  •  »b  ♦  m(0).  However,  the  more  precise  asymptotic  estimate 

(2.43)  yields  more  Information  concerning  continuity.  Using  (2.43)  and  formula  (9)  from 

Section  4,26  of  (31,  it  Is  not  difficult  to  show  that  the  solution  of  (1.1),  (1.2),  (1.3a) 

/ - - - 

is  continuous  across  the  lines  x  »  i  t  •  *b  +  m(0)  for  all  t  >  0,  and  the  solution  of 
(1.1),  (1.2),  (1.3b)  is  of  class  C1  across  these  lines  for  all  t  >  0. 
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E.  General  Remarks 


It  Is  clear  that  the  above  analysis  is  applicable  to  more  general  classes  of  singular 
kernels*  The  crucial  factor  is  the  asymptotic  behavior  of  m(X)  as  X  ♦  ",  which,  of 
course,  is  related  to  the  behavior  of  m(t)  as  t  ♦  0+.  If  m^(t)  -  mjtt)  as  t  *  0+ 

(and  some  mild  regularity  conditions  hold)  then  m^tX)  ~  m^tX)  as  X  ♦  •  in  any  sector  of 
the  form  |arg  X|  <  ♦  <  i.  (See,  for  example.  Theorem  33.2  of  18).)  Unfortunately,  such  a 
result  is  not  quite  sufficient  to  give  a  precise  regularity  theory  for  (1.1)  which  is  based 
solely  on  the  behavior  of  m(t)  as  t  ♦  0+.  However,  it  is  easy  to  see  that  our  arguments 
can  be  applied  to  any  kernel  whose  laplace  transform  is  suitably  behaved  in  the  right  half 
plane . 

If  m  belongs  to  L  (0,*)  and  la  completely  monotone  (lse#,  m  e  C  (<),«•)  and 
(-1)V«(t)  >  o  for  all  t  >  0,  k  «  0,  I,***),  then  m  can  be  continued  analytically  to 
the  slit  half  plane  Cn(-<“,0]  and  the  contour  of  Figure  1  can  be  used  in  (2.5)  for 
0  <  x  <  t  •k  +  m(0).  Moreover,  u,  u^,  u^,  ufct,  and  u^x  are  continuous  across  x  »  0 
for  t  >  0  and  the  jump  J  defined  by  (2.8)  is  analytic  on  (0,-).  The  only  nontrivial 
cases  in  which  J  vanishes  identically  occur  when  m  is  an  exponential  of  the  form 
(2.10).  (See  (2.13)  and  (2.15).)  We  conjecture  that  if  m  is  completely  monotone, 
integrable,  and  m(0+)  *  ",  then  the  solution  of  (1.1),  (1.2),  (1.3a)  is  continuous  (and 
the  solution  of  (1.1),  (1.2),  (1.3b)  is  of  class  C1)  across  the  lines 
x  «  ±  t  •  k  +  m(0 )  for  all  t  >  0. 

00  1 

It  is  to  be  expected  that  if  m  e  C  <0,»>  n  L  (0,“>,  then  the  solution  of  (1.1), 

„  / - - - 

(1.2),  (2.3)  is  of  class  C  away  from  the  lines  X  »  0,  x  t  •  *b  +  m(0).  However,  it 

does  not  seem  easy  to  prove.  A  result  of  this  nature  is  established  in  the  Appendix  under 

an  additional  assumption  on  m,  which,  rougly  speaking,  requires  that  derivatives  of  m 

are  not  "too  wild*  near  0. 
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3.  General  Initial  Value  Problems 

Using  the  results  of  the  preceding  section  and  the  method  of  superposition,  it  la 
straightforward  to  study  smoothness  properties  of  solutions  to  (1.1),  (1.2)  with  Initial 
conditions 

(3.1)  u(x,0)  "  u0(x),  u^ ( x, 0 )  ■  u^x),  x  e  R  , 

for  u0  and  u^  in  various  function  classes.  Indeed,  under  rather  mild  assumptions  on 
u0  and  uj,  the  solution  of  (1.1),  (1.2),  (3.1)  can  be  written  in  the  form 

(3.2)  u(x,t)  »  D*(x-5.t)u0(e)d£  ♦  {I*  uNx-C,t)u1(t)dt 

where  Ua  and  Ub  denote  the  solutions  of  (1.1),  (1.2),  (1.3a)  and  (1.1),  (1.2),  (1.3b), 
respectively.  (Recall  that  derivatives  are  to  be  interpreted  in  the  distributional  sense.) 

We  close  with  an  informal  discussion  of  the  regularity  of  solutions  of  (1.1),  (1.2), 
(3.1).  Our  assumptions  on  m  (and  b)  are  the  same  as  in  Section  2. 

A.  Regular  Kernels 

Suppose  that  m  satisfies  (2.6).  If  there  is  a  singularity  in  the  data  at  a  point 
x0  e  R,  then  three  waves  will  generally  emanate  from  this  point.  Singularities  will 
propagate  along  the  lines  x  “  Xg  ±  t  •  *  m(0),  and  (except  in  very  special 

circumstances)  there  will  be  a  stationary  singularity  at  x  «  Xg  for  t  >  0.  Such  a 
stationary  singularity  will  always  be  weaker  than  the  corresponding  propagating 
singularities. 

A  discontinuity  in  u^k^  produces  propagating  discontinuities  in  k»th  derivatives 

of  u,  while  a  discontinuity  in  u}  propagates  in  (k+1)-st  derivatives  of  u. 

ak+2u 

Discontinuities  in  Ugk*  or  u^'  both  lead  to  stationary  discontinuities  in  — 

3x 

If  m  is  a  single  exponential  of  the  form  (2.10),  there  will  be  cases  in  which  a 
singularity  in  the  data  does  not  lead  to  a  stationary  singularity  in  u.  See  Section  2A. 
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B.  Kernels  with  Power  Type  Singularities 


0* 

For  the  kernel  given  by  (2.17),  the  propagating  waves  will  always  be  of  class  C  . 

Consequently,  the  regularity  of  u  is  determined  by  the  stationary  singularities.  Discon- 

ak+2 

tinuities  in  ug  or  u}  both  lead  to  stationary  discontinuities  in  — — —  for 

3xlc+2 

t  >  0.  Thus,  there  is  a  definite  (but  finite)  gain  in  regularity. 

C.  Kernels  with  Logarithmic  Singularities 

The  situation  is  especially  interesting  for  the  kernel  given  by  (2.21).  Initially, 
the  propagating  singularities  are  stronger  than  the  corresponding  stationary  ones  and  the 
regularity  of  u  increases  in  a  manner  proportional  to  time.  This  will  continue  until  the 
strongest  stationary  singularity  and  the  strongest  propagating  singularity  are  of  the  same 
order.  Then,  the  stationary  singularity  will  dominate  and  the  regularity  of  u  will  not 

isiprove  any  further.  As  before,  discontinuities  in  Ug^ *  or  u^  lead  to  stationary 

,k+2 
o  u 

discontinuities  in  — rrr« 

3x* 1 

D.  Kernels  with  Weaker  Singularities 

For  the  kernel  (2.30)  (or  (2.41)),  the  propagating  singularities  are  stronger  than  the 
corresponding  stationary  ones  for  all  time.  The  smoothing  of  the  propagating  waves  leads 
to  only  an  infinitesimal  gain  in  regularity  of  u. 

For  the  kernel  (2.41),  we  can,  however,  make  a  definite  statement  concerning 
continuity.  If  ug  C  SV(tt)  and  u^  =  0,  the  solution  of  (1.1),  (1.2),  (3.1)  is 
continuous  on  R  x  { 0,») r  if  u0  I  0  and  ut  e  BV(R),  then  the  solution  of  (1.1),  (1.2), 
(3.1)  is  of  class  C1  on  R  x  (0,“) .  We  conjecture  that  this  is  also  true  for  the  kernel 
(2.30)  (and,  more  generally  for  any  completely  monotonic  m  e  L*(0,»)). 


Appendix 


In  this  appendix,  we  discuss  regularity  of  the  solution  u  of  (1.1),  (1.2),  (2.2) 


away  from  the  linee  x-0,  x  *  *  t  •  rb  ♦  m(0),  under  the  assumptions  that  m  is  non¬ 


negative,  nonincreasing,  and 
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SVetch  of  Proof:  The  argument  will  be  partitioned  into  several  steps.  Since  we  are 
interested  only  in  smoothness  (and  not  in  asymptotic  behavior  as  t  *  •>,  we  assume, 
without  loss  of  generality,  that  the  support  of  m  is  bounded. 

Step  t i  Consider  the  function  f  whose  Laplace  transform  is  given  by 
(a. 4)  f(X)  »  /b  -  m(l)  -  /B  . 

By  virtue  of  (a.1),  our  sign  conditions  on  m  and  b,  and  the  Riemann-Lebesgue  lemma,  it 
follows  from  Proposition  2.3  of  (16)  that  f  is  the  Laplace  transform  of  a  function 
f  e  L  (0,*).  He  want  to  show  that,  in  addition, 

(a. 5)  f  e  C^O,*)  . 

Observe  that  f  satisfies  the  equation 
(a. 6)  (f  +  /b  S)»(f  Wb  S)  »  b6  -  m  , 

or 

(a. 7)  (f*f)(t)  +  2/5  f(t)  »  -m(t)  , 

where  5  is  the  Dirac  distribution  and  the  *  denotes  convolution  on  [0,t).  We  can 
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where 


always  writs  m  in  tha  torn  n  ■  m,  +  m2, 

(a. 8)  m  (X)  “  0(— )  aa  1  *  «•,  Re  X  >  0 

1  Xn 

for  every  n  >  0  and  raj  nonnegative,  nonincreaaing,  and  satisfies 
(a. 9)  m2  e  L^O,-) 

with  Ira.l  .  arbitrarily  snail.  (Note  that  tmjtt)  e  t'co,®)  and  ltm'(t)l  »  In  I  .) 

L  1.  2  L 

OB 

It  follows  fron  (a. 7)  that  m1  e  C  [0,“),  and  consequently  derivatives  of  ra2  can 
grow  no  faster  than  derivatives  of  m  near  0.  It  also  follows  that 


(a. 10) 


/b  -  m(X)  -  /b  -  mj  (1 )  “  O(~0 


for  every  n  >  0  as  X  *  «*.  Therefore,  to  discuss  the  regularity  of  f,  it  suffices  to 
consider 

(a. 11)  f 2*f 2  +  2/b  f2  «  -n2  . 

By  choosing  ra2  sufficiently  small,  we  can  guarantee  that  (a.  11)  has  a  solution 

f2  e  1*^  (0,“)  with  Ifjl  1  arbitrarily  small. 

L 

We  rewrite  (a. 11)  as 

(a. 12)  2/1  f2(t)  +  2  /g/2f2<t«T)f2(T)dT  -  -  n2(t)  , 

from  which  we  find  formally 
(a. 13) 


2/b  If, I  „  -  2lf,l  „  •  If.l  <  Ira  I 

L  (t,-)  2  L  (t/2,«)  L  <0,->  L  [t ,- ) 


For  each  n  “  1,2,...,  let  us  set 


(a. 14) 


M  -  K-nln_l 

n  2  lVV) 


(a. 15) 


F  -  K~nl  f_l 
n  2  «•  -n  . 

L  [2  ,») 


where  X  is  a  suitable  positive  constant  (to  be  chosen  later).  It  follows  from  (a. 13) 
that 

08  OB  OO 

(a. 16)  2/B  T  F  -  2Klm,l  ,  I  F  S  T  M 

*■ .  n  2  ,1  *>  n+1  n 

n- 1  L  n“ 1  n»1 


If  K  is  large  enough,  then  \  M  <  * ,  and  (a. 16)  yields  a  bound  for 

n-1  n 

that  KImjl  y  is  sufficiently  small.  By  formally  differentiating  (a. 12) 


1  Fn  provided 
n- 1 

we  get  Jie  same 
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Kind  of  estimates  for  derivatives  of  itj  by  virtue  of  (a. 3). 

He  can  now  prove  that  fee  (0,")  as  follows.  Let  K  be  a  fixed  (but  arbitrary) 
positive  integer.  Hien,  we  write  m  «  m1  +  nvj,  with  re^  and  m2  as  above  and  mj  small 
enough  so  that  the  procedure  outlined  above  yields  estimates  for  derivatives  of  m2 
through  order  N.  For  each  e  >  0,  we  set 
(a. 17)  m2  (t)  -  uijlt+E),  t  >  0  , 

and  denote  by  f  the  solution  of 
e 

(a. 18)  2/5  f  (t)  +  2  j*/2  f  (t-t)f  <T)dT  =  -m  (t)  . 

e  £c  e 

It  is  clear  that  f2  6  C*  [(>,«•)  for  each  e  >  0.  Moreover,  our  estimates  for 

derivatives  of  m2  hold  for  derivatives  of  m2  -  independently  of  e.  As  e  ♦  0, 

e 

f2  ♦  f2  in  L^(0,*)e  in  addition,  f2  and  its  derivatives  through  order  N  are 
e  e 

1 

uniformly  bounded  on  compact  subsets  of  (0,*),  and  consequently  f 2  e  c  (0,“>).  Since 

a)  %r_  1 

f  -  fj  e  C  [0,«),  we  have  f  e  C  (0,*»).  Finally,  since  N  was  arbitrary,  we  conclude 
that  f  e  c“ (0,“), 

This  procedure  also  shows  that  derivatives  of  f2  satisfy  bounds  of  the  form  (a. 3), 

independently  of  the  particular  decomposition  of  m  into  m1  +  ra2-  Moreover, 

tf^(t)  e  L^O,")  with  ltf*(t)l  1  arbitrarily  small  if  ra2  is  sufficiently  small. 

L 

Step  2;  Next,  consider  the  function  g  whose  Laplace  transform  is  given  by 


(a. 19) 


g(\)  .  -L-  -  —  1  ■■ 

/b  -  all) 

1  1 
“  —  -  7 - -  . 

✓b  fMi  ♦  /b 


As  in  Step  1 ,  it  follows 
a  function  g  e  L^O,"). 
(a. 20) 
or 

(a. 21) 


from  Proposition  2.3  of  [16)  that 

We  note  that  g  satisfies 

-  gj*l/i  5  +  f)  -  S  , 
/b 


/b  g  +  f*g 


9 


is  the  Iaplace  transform  of 
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r 


Aa  before,  to  study  the  smoothness  properties  of  g,  it  is  enough  to  consider 
(a. 22)  Sb  q2  +  f2*g2  -  f2  • 

A  modified  version  of  the  argument  in  Step  1  shows  that  g  and  g2  have  the  same 

regularity  properties  as  f  and  f2<  respectively.  Moreover,  lg2l  1  and  ItgjitX  1 

L  L 

can  be  made  arbitrarily  small  by  choosing  m2  sufficiently  small. 

Step  3i  Let  c  be  a  positive  constant  and  consider  the  function  h  whose  Laplace 
transform  is  given  by 

(a. 23)  h(X)  ■  j  exp(cXg(X) ) 

2  * 

From  sign  conditions  on  the  kernel,  etc.,  we  know  that  h  e  L^qc [0,» ) .  Observe  that  h 
satisfies  the  differential  equation 

(a. 24)  —  h(X)  -  -  j  h(X)  +  (cXg(X))  •  h(X)  , 

which  translates  into 

(a. 25)  th(  t)  -  J*  h ( T ) dt  ♦  c  J*  Tg1  (T  )h(  t-T  ) dT  . 

As  in  Step  1,  to  study  the  smoothness  properties  of  h,  it  1b  enough  to  consider  the 
part  of  h  which  corresponds  to  m2,  i.e.  the  solution  h2  of 
(a. 26)  thj(t)  -  h2(T)dT  c  /*  xg' (T)h2(t-T)dt  . 

Indeed,  it  is  easy  to  see  that  h  -  h2  e  c”[0,«) .  The  argument  of  Step  1  cannot  be  applied 
directly  to  (a. 26)  because  of  the  singular  coefficient  on  the  left-hand  side.  To  overcome 
this  difficulty,  we  fix  tg  >  0  and  set 

(a. 27)  (t)  -  h2 (t+tg ) ,  t  »  0  . 


After  several  changes  of  variable,  (a. 26)  yields 


(a. 28)  <t+t  )*  (t)  “  /o  (T)dT 

0  0 

+  c  ll/2  Tg^(T)*t  (t-T)dt 
0 

♦  C  f*/2  Tg*(T)*t^(t-T)dT 
t 

+  /0  h2<TXJT 
t- 

+  c  /0°  (t-t-t0-T)g'(t+t0-T)h2(T)dT  . 

We  know  a  priori  that  $t  ,  h2  e  L^tO,*).  Moreover,  we  have  information  on  g2  fro“ 
Step  2.  A  modified  version  of  the  procedure  of  Step  1  can  now  be  used  to  show  that 
h  e  c”(0,«>).  (In  this  case  intervals  of  the  form  [2“n,T)  with  T  <  •  should  be  used  in 
place  of  [2  ”,“).) 

step  4 ;  For  a  -  1,  8  «  0,  the  integrand  in  (2.5)  can  be  written  in  the  form 
(a. 29)  exp[l(t  -  ^]J  •  j  •  exp(xXg(X))  . 

(The  case  of  general  a,  8  follows  easily  from  this  one.)  We  can  therefore  conclude  from 
step  3  that  for  each  (fixed)  x  >  0,  u  is  a  c“  function  of  t  -  x//B  for 
t  -  X/Vn  >  0,  and  consequently  u(x,*)  e  C  (x//b,“). 
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